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ABSTRACT 

Recent claims that the strength I?igmf of the intergalactic magnetic field 
(IGMF) is > 10 -15 G are based on upper limits to the expected cascade flux in 
the GeV band produced by blazar TeV photons absorbed by the extragalactic 
background light. This limit depends on an assumption that the mean blazar 
TeV flux remains constant on timescales > 2(I? I gmf/10~ 18 G) 2 /(E/10 GeV) 2 yr 
for an IGMF coherence length « 1 Mpc, where E is the measured photon energy. 
Restricting TeV activity of 1ES 0229+200 to ~ 3 - 4 years during which the 
source has been observed leads to a more robust lower limit of Bjgmf <i 10~ 18 
G, which can be larger by an order of magnitude if the intrinsic source flux above 
« 5 - 10 TeV from 1ES 0229+200 is strong. 

Subject headings: gamma rays: theory — radiation mechanisms: nonthermal 



1. Introduction 

The measurement of the intergalactic magnetic field (IGMF) gives information about 
processes operating in the early universe that are imprinted on the large scale structure of 
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the universe (see iNeronov &: Semikoa 120091 ) . Faraday rotation measurements of the radio 
emission of quasars, patterns in the arrival directions of UHEC Rs towards the supe rgalactic 
plane and Cen A, and theoretical arguments from COBE data (IBarrow et al.lll997l ) indicate 
that -Bigmf ^ 10~ 8 G, but no direct measurements or lower limits of the IGMF in the 
voids have been firmly established. Gamma-ray a stronomy pr ovides a method to measure 
the IGMF through magnetic field-induced delays (lPlagalll995l). pair ha los from sources of 
TeV photons direc ted away from our line of sight ( Aharonian et all 1 1994), and halos around 
jElviv et al.ll2009l) and cascade-radiation spectra of (iDai et al.l I2OO2I ; iD'Avezac et al.l I2OO7I : 



Murase et al 



20081 : INeronov fc Vovkll2010l : iTavecchio et al 



2010al lbh point sources of high- 



energy 7 rays. 

In the latter approach, TeV photons from cosmic 7-ray sources interact with photons 
of the extragalactic background light (EBL) to produce e~e + pairs. If the pairs are not 
significantly deflected by the IGMF, cosmic microwave background (CMB) photons, which 
dominate the radiation energy density in intergalactic space, are Compton-scattered in the 
original direction of the pairs. Deflection out of the TeV beam depends on the jet opening 
angle 9 y By comparing blazar TeV fluxes with upper limits on the GeV radiation flux 
measurements with the Fermi Gamma ray Space Telescope, sever al claims (see Table 1) 



have been made that lo w er limits on the IG MF have been measured (INeronov fc Vovkl 12010 



Tavecchio etaD boiOal Jbl: | Dolag et al.l 120 111 ). These limits, which depend on the assumed 
opening angle 9- } of the TeV photon source, are summarized in Table 1 for the TeV blazar 
1ES 0229+200 at redshift z = 0.1396. Under the assumption of persistent TeV emission 
over long time scales, these studies find that -Bigmf ~ 10~ 15 G for magnetic coherence (or 
correlation) lengths \ co h ~ 1 Mpc when 0- s = 0.1. 

Implicit in all these studies is that the TeV blazars used to infer the IGMF emit constant 
flux over a long period of time. Because blazars are highly variable, a more defensible limit 
is obtained (lacking other ways to infer a blazar's TeV activity lifetime) by assuming that the 
TeV emission is emitted only over the past few years during which it has been monitored. 
A simple semi-an alytic approach is used to derive min imum values, for comparison with 
numerical models ( IDolag et al.l 1201 It iTaylor et al.l 120 111 ). 



2. Time Delay and Deflection Angle of Emission from First Generation Pairs 

Consider a source and observer separated by distance d, as shown in Figure 1. Photons 
with dimensionless energy e\ = hv\jm e (? = 2 x 10 6 i?i(TeV) emitted at angle Q\ to the line 
of sight between the source and observer, travel a mean distance A 77 = A 77 (ei,2:) before 
converting into an electron-positron pair via 77 absorption with photons of the EBL. The 
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Fig. 1. — Sketch of the geometry of the process. A photon with energy -EreV TeV, emitted at 
angle 9i < 9 S to the line of sight, interacts with an EBL photon to create an electron-positron 
pair with Lorentz factor 7 = 10 6 7 6 . The lepton is deflected through angle # dfl and scatters 
a CMB photon to energy -EgcV GeV, which is observed as a source photon by the Fermi 
LAT if it is detected at an angle 9 < 9 ps f(EQ e y) to the source. The underlying simplifying 
kinematic relation in the semi-analytic model is 76 ~ -Erev ~ V^GeV- 



pairs cool by scattering CMB photons to -EceV GeV energies, which are detected at an angle 
9 to the line of sight to the source when the secondary electrons and positrons (hereafter 
referred to as electrons) are deflected by an angle # dfl . The GeV emission, to be detected, 
must be within the energy-dependent Fermi LAT point-spread-function an gle 9 ps f(EQ e y). 



The system is treated in the low redshift limit (cf. iNeronov Sz Semikozll2009l ). 



The time delay At between the reception of photons directed towards the observer and 
those formed by the process described above is given by 



<isin( 



cAt = A 77 + x — d = A 77 H — d = 

sin 9 d Q 

A 77 (l — cos 6^) — d{l — cos 9) , (1) 

noting that x = dsm9\/ sin^fl and A 77 = dsin.9/ sinfj'dfl. In the limit of small observing and 
deflection angles, equation (0Q) implies 

At = ^ figa , (2) 
provided that the photon is detected at an angle 

9 = A ^ (ETcv) ^ /KEGev) < 9 psf (E GeV ) (3) 

to the source. Note that the deflection angle depends on either the primary photon energy 
E TeV or Compton-scattered photon energy E GeV , since they are related by E Ge v ~ E^y, as 
we now show. 

The average CMB photon energy at low redshift is €q w 1.24 x 10 -9 in m e c 2 units, so that 
mean Thomson-scattered photon energy is ex ~ (4/3)e 7 2 , where 7 = E Te y / \2m e c 2 ) implies 
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76 = 0.98i?TeV- Thus, an electron with Lorentz factor 7 scatters CMB radiation to photon 
energy E when 7 6 = E Te y — l-lV^GeV- The characteristic length scale for energy losses due 
to Thomson scattering is At = 3m e c 2 / '4<7t«cmb7 = (0.75/7 6 ) Mpc, where uqmb — 4 x 10~ 13 
erg cm -3 is the CMB energy density at low redshifts. While losing energy, the electron is 
deflected by an angle 9b — At/Vl in a uniform magnetic field of strength I?igmf = 10~ 15 -B-i5 
G oriented perpendicular to the direction of motion of the electron, where the Larmor radius 
r L = m e c 2 'j /eB = 0.55(7 6 /-B_i 5 ) Mpc. Thus the deflection angle for an electron losing energy 
by scattering CMB photons to energy E in a uniform field is 9b = At/?"l — I.IB^/Eq^. 
Introducing a coherence length X co h that characterizes the typical distance over which the 
magnetic field direction changes by ~ tt/2, then the deflection angle 



if At < X C oh 

9 da ee w9 B , with w = { /a^ if At > XcQh _ . (4) 



For 1ES 0229+200, TeV radiation has been detected to energies E < 12 TeV ( lAharonian et al 



20071 ). with an « 15% error in the energy measurement. An uncertainty in the analytic treat- 



ment is that the mean free path A 77 (-Erev) varies by a factor of ~ 2 between z — > and 



z = .14, and between different EBL models. For instance, the EBL model of iFinke et al. 
fcoioh gives A 77 (E) ^ 200 Mpc, 125 Mpc, and 7 Mpc at E = 1, 3, and 10 TeV, respectively 



and a low EBL model based on galaxy counts (IKneiske &: Dold 120101 ) gives A 77 (£ l ) = 280 
Mpc, 150 Mpc, and 85 Mpc, respectively. For analytic estimates, we write A 77 = IOOA100 
Mpc, though we use the accurate energy dependence of A 77 (-&r e v) i n the numerical calcula- 
tions. The importance of pair-cascade radiation with angular extent broader than the Fermi 
LAT psf depends on the value of 

V/ ^ d9 ps f(E Gc y)/9 df i ^ T 11 {E Tc y)9 ps f(E GeV ) . . 

A 77 A 77 9^i 

where A ps / is the effective distance a primary photon would have to travel to make a GeV 
photon detected at the edge of the Fermi psf given the parameters of the IGM. The value 



Table 1: Derived Limits on B IGMF for the source 1ES 0229+200 



1ES 0229+200 


4 (rad) 


-Sigmf(G) 


Neronov & Vovk (2010) 
Tavecchio et al. f2010a) 
Tavecchio et al. (2010b) 
Dolae et al. (2011) 


71 

0.1 

0.03 

0.1 


> 3 x 10" 16 

> 5 x 10~ 15 

> 2 x 10~ 15 

> 5 x 10~ 15 
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of 9 ps f(EG e y), taken here as the 95% Fer mi L AT confinement angle, is f rom the Fermi 
LAT instrument performance paga__ (see also Rando 2009 ; Taylor et al. 2011 ). For the EBL 
model of Finke et al. (2010), the cascade emission can be treated as a point source when 
B_ l5 < 0.05£° e 6 v for 0.2 < E GcY < 20. 

For a source at distance d = dc P c Gpc, with dc pc ~ 1 corresponding to z ~ 0.2, the 
time delay for emission observed at angle 



6 = 0.01 



A 



100 



B-i 5 w 



d Gpc K E/10 GeV y 
from the line of sight is given from equation (j2j) by 

B-i 5 w 2 



At(yr) = 2 x 10 A 



100 \ 



E/10 GeV y 



(6) 



(7) 



Short delay times are restricted to conditions of small -Bigmf and large E where, as just 
seen, extended pair halo emission can be neglected. 



Equation (j7|) shows that small time delays are implied when A 77 is small and X ps f/\ 



77 



> 



1. When A 77 < At, an additional delay ~ X^O^/c ar ises during the time that the electrons 



are losing energy and being deflected by the IGMF ( IMurase et al.l l2008t llchiki et al.l 12008 



Razzaque et al.l I2004J ) . Such small values of A 77 ~ 1 Mpc are only relevant at low redshifts 
for > 100 TeV photons that pair-produce within w 1 Mpc of their source, where the magnetic 
field may not be representative of the dominant volume of the voids. 



3. 7 Ray Data of 1ES 0229+200 



The TeV blazar 1ES 0229+200, which provides some of the st rongest constraints on the 
lower limit to the IGMF field strengt h, was observed with HESS (lAharonian et al.l 120071 ) in 
2005 and 2006 and with VERITAS JPerkins I bold ) in October 2009 - January 2010. No 
evidence for variability of the TeV flux has been reported, so the observations give an average 
TeV flux from this source on timescales of ~3 yr, though with poor sampling. The HESS 



and preliminary VERITAS data ( jPerkins 1120101 ) are shown in Fig. 2 by the blue open circles 
and red squares, respectively. 

Fermi Large Area Telescope upper limits on TeV blazars were reported previously 
flAbdo et al.l bood . boioh . Here we reanalyze the Fermi LAT data for 1ES 0229+200 col- 
lected from 2008 August 4 to 2010 September 5 in survey mode. To minimize systematics, 



www-glast.slac.stanford.edu/software/IS/glastJat_performance.htm 
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only photons with energies greater than 100 MeV were considered in this analysis. In or- 
der to avoid contami nation from Earth-l imb 7 rays, a selection on events with zenith angle 
< 105° was applied ( lAtwood et al.ll2009l ). This analysis was performed using the standard 
likelihood analysis tools that are part of the Fermi Science Tools software package (version 
v9rl5p5)@ The P6_V3_DIFFUSE set of instrument response functions was used. Photons 
were selected in a circular region of interest (ROI) 10° in radius, centered at the position 
of 1ES 0229+200. The isotropic background, including the sum of residual instrumental 
background and extragalactic diffuse 7-ray background, was modeled by fitting this com- 
ponent at high galactic latitude (isotropicJem_v02.txt, available from the FSSC website). 
The Galactic diffuse emiss ion model version "gllJem_v02.fit," was used in the analysis. The 
profile likelihood method ( IRolke et al.l 120051 ) was used to extract 95% confidence level upper 
limits at the location of 1ES0229+200 assuming a power-law energy distribution with pho- 
ton index=2, all 1FGL point sources lying within the ROI being modeled with power-law 
distributions. The upper limits shown in Figure 2 are obtained in the energy bins 0.1-1 
GeV, 1 - 3 GeV, 3-10 GeV, 1 - 10 GeV, and 10 - 100 GeV. 



Model for Cascade Radiation 



The limits on the IGMF can be established by employing a simple semi-analytic model 
for the cascade radiation spectrum. Using the notation that f t = vF v at dimensionless 
photon energy e, and that each photon is attenuated into a pair with each electron taking 
one-half the original photon's energy, then a straightforward derivation gives 



-(-) 2 
2 V 



A\/ es/4e ,7dfl,7(A*en fl )] 



g?7 7" 



47V 



IX 



/ £ {exp[r 77 (e,z)] - 1} 



where 7^ = e/2. The interior integrand represents the fraction of deabsorbed source photon 
flux converted to pairs, an d the exterior integral represents the Compton-scattered spe ctrum 
from cooled electrons (cf. iRazzaque et al.ll2004l : iMurase et al.ll2008l ; llchiki et al.l 120081 1 . The 
opacity due to EBL attenuation for photons with measured dimensionless energy e from a 
source at redshift z <C 1 is T 77 (e, z), and depends on the EBL model. 



Equation (JSJ) employs the isotropic Thomson kernel, with the CMB radiation approx- 
imated as a monochromatic radiation field, but the results in Figure 2 are also integrated 



2 http://fermi. gsfc.nasa.gov/ssc/. 
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over the energy distribution of the blackbody radiation field. The use of the KN kernel makes 
negligible difference for photons with energy < 20 TeV. In the three terms in the lower limit 
of the exterior integration, the first gives the kinematic minimum electron Lorentz factor to 
scatter a CMB photon to energy e s . The second is the value of the deflection Lorentz factor 
7dA obtained by equating the Thomson cooling time and the timescale Qfi^jc when the elec- 
tron is deflected outside the photon beam of opening angle 6 y The third limit, 7(At eri9 ), 
represents the Lorentz factor to which electrons have cooled after the blazar engine has 
been operating for time At eng , and follows from equation fl2]) by solving At(7 en5 ) < At eng for 
leng = l(At eng ). Here we approximate A 77 (i?Tev) ~ d/T 77 (E^ e y) Mpc, using a fit to the Finke 
et al. (2010) EBL model for 1ES 0229+200. A calculation with A 77 (£ TeV ) « d/(2r 77 (£ TeV )) 
Mpc gives similar results. Only the first generation of cascade emission attenuated by the 
factor exp[— r 77 (ei, z)] is shown here. 

Results of calculations using this simplified analytic model are shown in Figure 2. Fig. 
2a is a calculation where the blazar engine operates for indefinitely long times, with the 
reduction of cascade flux due to deflection away from the beam for a jet and the detection 
of a plateau flux of iso tropized radiation determined by the jet opening angle 6*j = 0.1 
( ITavecchio et al.l l2010af ) . The source spectrum is described by a super-exponential cutoff 
power law vF v oc £ 4/5 exp[— (B/5 TeV) 2 ] in Figs. 2(a) and 2(b), and by an exponential cutoff 



power law vF v oc E 4 / 5 e xp(— E/10 TeV 



in Fig. 2(c). In agreemen t with previous results 



( INeronov &: Vovkl l2010t ITavecchio et al. 

3 x 10 -16 G is implied in order to reduce the GeV flux below the Fermi upper limit 



2010aJJbJ; iDolag et al.l 1201 if ), a value of B 1GMF > 

From 



the calculations, we also find that under the assumption of persistent TeV blazar emission, 
halo emission becomes increa singly dominant fo r large jet opening; angles. Dete ction of halos 
around AGNs, as claimed by lAndo fc Kusenkol (120101 ) (cf. INeronov et al.l 1201 it ), would then 
favor detection in sources with large opening angle, long lived TeV engines. Also under the 
persistent emission hypothesis, a maximum jet opening angle 9- } < 0.4 is implied in order 
that the isotropized radiation does not violate the Fermi LAT upper limits. 

The effects of -Bigmf on the received spectrum of reprocessed TeV radiation when the 
blazar engine is assumed to emit a constant TeV flux over an engine time At eng = 3 yr 
are shown in Figs. 2a and 2b. These calculations show that -Bigmf ^ 3 x 10~ 19 G for the 
case where the assumed source spectrum is sharply cut off above 5 TeV. Uncertainties in 
the analytic model, including the strong sensitivity of the cascade spectrum on 7 en9 , relaxes 
our conclusions to an analytic, order-of-magnitude minimum IGMF of -Bigmf ^ 10~ 18 G for 
At eng = 3 yr. Fig. 2c shows that the minimum magnetic field also depends sensitively on 
the characterization of the high-energy spectral flux, which can then q uickly cascade in to 
the 10 - 100 GeV band and violate a Fermi upper limit (or detectio n; seelOrr et al.l (120111) ). 
By assuming source spectra with larger fluxes above « 5 - 10 TeV, IDolag et al.l (1201 if ) and 
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Taylor et al.l ( 1201 ll ) derive larger values for the minimum -Bigmf , but not more than a factor 



of a few above the analytic results when difference in activity times and primary source 
fluxes are considered. 



5. Discussion and Summary 



Previous GeV/TeV inferences of the strength of the intergalactic magnetic field make 
an assumption that the mean blazar TeV flux over millions of years remains similar to values 
observed over the last few few years. Our knowledge of the blazar engine is not yet so good as 
to have high confidence in this assumption, though some models for slowly varying TeV flux 
from TeV blazars can be noted. For example, a slow cooling rate of the electrons that make 
the TeV photons could imply a slowly varying 7-ray flux even if the blazar engine is very 
active. For electrons scattering photons to TeV energies, the synchrotron cooling time for the 
obse rver is t Rim = ( 1 + z)Q 7Tm P c/ (SnCrB' 2 ^) = 50/i?(TeV) yr, using the fitting parameters 



of iTavecchio et al.l (l2010d ) for 1ES 0229+200 (break Lorentz factor 7^ = 5 x 10 6 , emission 



region magnetic field B' — 5 x 10 4 G, and Doppler factor <5d = 40). Relativistic electrons in 
an extended jet that Compton scatter photons of the CMB c ould also make slowly varying 
TeV radiation in sources like 1ES 0229+200 or 1ES 1101-232 faottcher et al.ll2008h . In this 



model, relativistic electrons lose energy on timescales of « 750/[(r/10) V^TW)] yr . These 
models do not, however, provide good reasons to expect TeV blazars to produce steady flux 
for thousands or millions of years. 

A more reliable limit is obtained from direct measurements of TeV fluxes. For the 



hand f ul of observati ons of 1ES 0229+200 over 3-4 years of observing ( lAharonian et al. 
20071 ; iPerkins 1 12010| ). no TeV flux variations have been reported. Using such timescales 
leads to a limit of Bi GMF (G) > 10~ 18 (£'/10 GeV)y / At/3 yr/\/X W o, assuming that X coh « 1 
Mpc. By assuming strong intrinsic > 10 TeV emission from 1ES 0229+200 (which is not 
observed because of EBL attenuation), Fermi LAT flux upper limits at ~ 100 GeV can 
be violated, leading to larger limiting values of -Bigmf (G) > 5 x 10~ 18 G. Evidence for a 
strong primary flux at > 10 TeV co mes from detection of a shoulder feature at ~ 1 TeV, as 
found in the numerical calculations (jDolag et al.l 1201 if ) and analytical results (Fig. 2c), and 
suggested by the joint VERITAS/HESS data. Note that our calculations assume negligible 
contribution from cascades induced by photopair interactions by > 10 18 eV cosmic rays 
( lEssey et al.ll2010l ). More frequent, sensitive, and broadband GeV - TeV observations of 
1ES 0229+200 can test whether the average TeV flux corresponds to the flux that has been 
historically measured or is unusual. 



Evidence for long-lived TeV radiation can be identified in pair halos ( lAharonian et al 
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1994 ) from misaligned blazar candidates such as Cen A or M87. Searches for pai r echoes 
from GRBs, which are sensitive at <C lCT 21 /A co h(Mpc) G (ITakahashi et al.l 120081 ). would 
test our claim that -Bigmf ^ 10~ 18 G. A lar ge field-of-view detector like the High Altitude 
Water Cherenkov telescope (IGoodman Il2010l ). or systematic monitoring campaigns of blazars 
like 1ES 0229+200, 1ES 1101-232 [z = 0.186), 1ES 0347-121 (z = 0.185) or other bright, 
moderate redshift BL Lacs with the present generation of air Cherenkov telescopes or an 
advanced Cherenkov telescope array, will give better information about the duty cycle of 
TeV blazars and provide more secure constraints on the value of the intergalactic magnetic 
field. 
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Fig. 2. — Model of cascade radiation spectrum, equation (jSJ), applied to HESS, VERITAS, 
an d Fermi obse r vation s of 1ES 0229+200, using model spectra (solid curves) and EBL model 
of iFinke et al.l ( 120101 ) to give attenuated source spectrum (dotted curves), (a) Cascade 
spectra for 1ES 0229+200 assuming persistent TeV emission at the level observed with HESS 
and VERITAS, for different values of -Bigmf and = 1 Mpc (solid) or X coh = 100 kp (dot- 
dashed). The psf constraint for the \ co h = 1 Mpc case is shown by the dashed curves. 
Cascade spectra when source radiates TeV flux for 3 yr with constant average spectrum 
given by power-law with vF v index = 4/5 for source spectrum with superexponential cutoff 
oc exp[— (E/5 TeV) 2 ] (b) and exponential cutoff oc exp(— (£7/10 TeV) (c), are shown for the 
case X co h = 1 Mpc with different values of £>igmf, as labeled. 



